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If the spontaneous breaking of Peccei-Quinn symmetry comes from soft supersymmetry breaking, 
the fermionic partners of the symmetry-breaking fields have mass of order the gravitino mass, and 
are called flatinos. The lightest flatino, called here the flaxino, is a CDM candidate if it is the 
lightest supersymmetric particle. We here explore flaxino dark matter assuming that the lightest 
ordinary supersymmetric particle is the stau, with gravity-mediated supersymmetry breaking. The 

. . . decay of the stau to the flaxino is fast enough not to spoil the standard predictions of Big Bang 

OO ' Nucleosynthesis, and its track and decay can be seen in future colliders. 

o , 

O ! PACS numbers: 11.30.Hv, 12.15.Ff, 12.60.Jv, 14.60.Pq, 14.80.Mz, 98.80.Cq 

I. INTRODUCTION 

■ By observing the temperature anisotropy of Cosmic Microwave Background (CMB), the WMAP Collaboration 
' reported that the density parameter of dark matter (DM) or cold dark matter (CDM) riDM,obs is 



Oh' 
I 



o . 



o 

00 

o 



^^M,ohsh' = 0.104t°:J]°^, (1) 



at 68 % C.L. [H where h is the normalized Hubble constant. One of the most popular candidates for DM is the 
CIh' lightest supersymmetric particle (LSP) in supersymmetric theories with R-parity. Recent studies showed that Big 
1^ , Bang Nucleosynthesis (BBN) puts interesting limitations on the candidates of LSP DM. 

A long-lived next LSP (NLSP) whose lifetime is longer than 0.1 sec is dangerous for the successful BBN because the 
] decaying NLSP produces a lot of (high-energy) daughter particles such as photons and hadrons during/after BBN, 
CNJ ■ which can destroy ''He (and D) and non-thermally produce another light elements copiously 0, d, Q- In particular, 
^ ' a non-thermal hadron emission due to the decaying NLSP is severely constrained by the observational light-element 
25 ■ abundances of D, ^He, ^Li, and ^Li [i,i3. 

In the popular case of the neutralino (gaugino or higgsino) LSP, as is well-known, the unstable gravitino decay 
puts stringent upper bound on the reheat temperature after primordial inflation, which could be problematic for 
various cosmological considerations, e.g., of the inflation models. This problem can be avoided if the gravitino is the 
LSP. (See also [5| for constraints on the reheat temperature in the gravitino LSP scenario). However, the decays of 
ordinary supersymmetric particles to the gravitino turn out to be dangerous as well. The case of neutralino NLSP, 
because the branching ratio into hadrons is close to the order of 0(1), is almost excluded [1, 0| by considering BBN 
• • constraints [i,[3. SleptonNLSP scenarios appear much more favorable than neutralino NLSP scenarios in cosmology. 
Although the sneutrino NLSP scenario should be possible because of the smaller hadronic and radiative decaying 
modes, the primordial abundance of the sneutrino is too small to satisfy the dark matter density for its daughter 
j_j LSP particle I, i]. Therefore, we would need another production mechanism of the LSP in cosmology, which usually 
5^ , requires a fine tuning. 

Recently, the stau NLSP scenario has been discussed a lot as there are specifically severer/more interesting problems 
in BBN related to the long-lived negatively charged particle [l^, [HI [12, [ll| • See also long-term problems discussed in 
literature [14J and topics in doubly-charged particles • If the lifetime of the negatively charged particle C~ such as 
f~ is longer than 10^ sec, bound states with ambient light elements can be produced. In particular, the bound state 
with "^He denoted by ("^HCjC"), which is produced after 10^ sec, significantly enhances the reaction rate D-|-('*He,C~) 
^^U + C- +-f^hya factor of ~ 0(10^) [11]. Then ^Li is over-produced by the catalyzed BBN, and this scenario 
is severely constrained by observations [l7| . 

Various attempts [Tsl [Tol . [20l . [2]| have been made to tackle this problem and realize a successful dark matter 
production mechanism. In this paper, we suggest the "flaxino" DM as a viable option to avoid the above-mentioned 
difficulties for the neutralino or slepton NLSP. The flaxino is an axino-like fermion appearing in flat-direction axion 
models where the spontaneous breaking of Peccei-Quinn symmetry comes from soft supersymmetry breaking and thus 
all of the fermionic partners of the symmetry-breaking fields have mass of order the gravitino mass. The flaxino has 
a coupling to the ordinary superparticles suppressed by the axion scale Fa and thus attractively leads to a fast decay 
of NLSP, which should be compared with the gravitino coupling suppressed by the Planck scale. 
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II. LSP FROM FLAT-DIRECTION AXION SECTOR 



The strong CP problem can be naturally resolved by assuming the Peccei-Quinn U{1) symmetry spontaneously 
broken at the usual high scale, Fa = 10^°'^^ GeV, and thus leads to an axion of the KSVZ [2j| or DFSZ type [H. In 
the supersymmetric standard model, the realization of the DFSZ axion model nicely leads to the origin of the term 
[23 |. Furthermore, the axion scale Fa appears as the geometric mean of the Planck scale and soft supersyrnmetry 
breaking scale by invoking an almost flat potential for the Peccei-Quinn symmetry breaking fiel ds [251 [2a |. For 
instance, let us introduce two singlet fields P, Q with appropriate Peccei-Quinn charges to allow [l^, [23 

W^,^^H,H, + ,gi (2) 

where Mp is the reduced Planck scale. Then, the /i term and axion scale can be related by /i = h{P){Q)/Mp and 
Fa ~ (P), (Q) y/rfiMp - 10^° where m ~ lO^-^ GeV is the typical soft mass scale (correctly Fa = ^(P)^ +9{Q)^ 
in this model). This is below the scale Fa ^ 10^^ GeV required to make the axion the CDM. This model contains two 
scalar flatons Pi_2jOne pseudoscalar fiaton F' , and two flatinos Pi,2- This kind of flat DFSZ axion models lead to 
thermal inflation 29] and the model parameters are rather severely constrained in order to avoid the over-production 
of unwanted relics like axions or the LSP after thermal inflation ^2^. An interesting collider implication of such 
models has been pointed out in Ref. 30]. The LSP can be some linear combination of the fermionic parts of P and 
Q, which we now call "flaxino" (denoted by Fi), and then the usual neutralino or stau NLSP decays to the flaxino 
mainly through its mixing with higgsinos driven by the ^ term. Although the corresponding decay lengths can be 
larger than the size of detectors, the copious production of the NLSPs in future colliders may enable us to observe 
their decays for the low axion scale Fa ^ 10^° GeV [s^l- The general idea of the axino DM has been put forward in 
Ref. 31] (see also a recent related paper ^32]). 

Some of the current authors have investigated the scenario of the "flaxino" LSP and the neutralino NLSP in the 
flat DFSZ axion models using gravity- mediated SUSY breaking [13, [11] . The axino LSP in KSVZ type models has 
been considered in Ref. [33] . The KSVZ type models (hadronic axion models) do not agree with the framework of 
the flat-direction models because of the overproduction of the flaxino abundance, which is excluded by BBN [27j . 
Extending the analysis of Ref. [1^ where the neutralino NLSP is assumed, we will consider the possibility of the stau 
NLSP which is also allowed in a large parameter space of the MSSM with the minimal supergravity scheme. 

Let us first note that the DFSZ models predict a shorter lifetime of the stau NLSP compared to the KSVZ models 
since the stau-axino/flaxino coupling arises at tree level in the former case and at one- loop level in the latter. As we 
will see later, even in the case of Fa ^ 10^"* GeV, the stau lifetime does not exceed 10^ sec for typical mass scales of 
SUSY particles in the DFSZ axion models and does not induce the ^Li over-production problem. This is a distinct 
feature of the DFSZ axion models. Recall that the usual upper bound, Fa ^ 10^^ GeV, can be relaxed in late-time 
entropy production scenarios such as thermal inflation models which we are considering at present. 

From the viewpoint of detectabilities of a superparticle at Large Hadron Collider (LHC), such a long-lived charged 
particle is promising as it can be easily traced within the detector. Moreover, quite recently Refs. [34, 35, 36] have 
discussed aggressive future plans to place new stopper, e.g., 5-10 meter wall made of iron, water tanks, rocks and 
so on, outside regular detectors in LHC such as ATLAS or CMS to electromagnetically stop the produced staus. If 
these plans are realized, we will be able to measure the stau's lifetime even if its decay length is longer than the size of 
the detectors (^ 0(10) m), and reconstruct the track of the neutral LSP by catching other daughter particles. Then, 
according to a similar discussion in Brandenburg et al's analysis [3^, observing three-body decay of the staus will 
enable us to distinguish between signals of the axino/flaxino LSP and the gravitino LSP. 

We refer to [13] for another topics of axino DM and gauge-mediation supersyrnmetry breaking models, and Ref. [11] 
for general cosmological constraints on a decaying saxion or flaton. See also [39j] for an inverted mass case such as 
the axino NLSP and the neutralino LSP, and [4^] for the gravitino NLSP and the axino LSP in gauge-mediated 
supersyrnmetry breaking models. 



III. STAU LIFETIME 

Following Ref. [s^], we calculate the lifetime of the stau NLSP which decays into the flaxino LSP Fi in the DFSZ 
axion models. For the sake of simplicity we show only mixing between Fi and the usual four neutralino components; 
bino, wino, and two higgsinos. This case is sufficient if we are not interested in the heavier flaxino F2 whose mass is 
larger than m ^0 . Then, we get 

r(n - rFi) = ^v'A(l,r2,r2)[(l4^ \' + \bl^ \') (l - - r^) ~ 4ror.Re (4%^*)] (3) 
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where ruf-^ is the hghtest stau mass (= rrif in this paper), tq = m^^/mfj, = 7X1^-1 m^^ with 

A(a::, y, z) — + + — 2xy — 2yz — 2zx, (4) 

and 

flo' = V2sj: [ggT3fNo2 + g {q/ - TsfNg^)] - c^gN^^nir/ V^cpmw, (5) 

= V2cj:g'qfNoi + SfgNt^s.mr / V^cpmw ■ (6) 

Here m^. and tovv' ^re the masses of the tau lepton and the weak boson, g and g' are the weak and electromagnetic 
couphng constants, respectively, and the Z-boson couplings to quarks and leptons for the tau lepton are give by 
(T^^^qr) = (—1/2,-1). The matrix elements Noi (i=l,2.3,4) are flaxino-neutralino components of a unitary matrix 
Nij which is the 5x5 mass matrix of four neutralinos plus one flaxino, M^f (fc, Z=0, 1,2.3,4) in the {Fi,B, W°, i?^, H^) 
basis. It can be diagonalized as 

N;kMi?Ni,^S.,,mj^. (7) 



where 



/ mp^ -Sispn -5iC0ii \ 

Ml —ci3Swmz spswrnz 

M2 cpcwTUz —SjscwTnz 

—SiSjSiJ, —CfjSwmz CfjCwmz — /i 

\ -SiCf3H SfjSwmz -spcwmz / 



(8) 



Here sp = sin/3, Cf3 = cos/3, sw = sinOw, and cw = cosOw with the Weinberg angle dw, and 



61 = — '^^ f^Qg (f) + X sin 4>] , (9) 

Fa X \ J 

with the electroweak scale v — 264 GeV, x = {P)/{Q), the mixing angle cj) between Fi and F2 which is expressed by 
using X to be cos 20 — —Xj^Jx^ + 1 and sin 20 = —x/\/x'^ + 1 [l^. The mass eigenstate labels follow the convention: 
mp_^ '^rrifj^^ < < rrifj^ < m^^ < m^^ (Note that m p^ can be larger than m although we did not write it down 
explicitly in ([5])). In the stau sector, the diagonalization is taken as follows; 



(10) 



Because the mixing angle in and Cf is an unknown parameter, we will study the Sf dependence by changing it 
from zero to unity. 

Note that the stau decay rate is approximately represented by 



r(fi^rFi) ~ ^ 



di- 
ll 



(Im X A) 



2 1 

X ^ (11) 

-1 



for Fa ^ 10^" GeV and rrif ^ jj. ^ 0(10^) GeV with a suppression factor A ~ 0(10) caused by electroweak couplings, 
mixing angles, kinematic suppressions and so on. Thus, we can expect the decay signals of staus within the meter-size 
detectors in LHC. 

IV. ABUNDANCE OF THE FLAXINO DM 

The flaxino LSP abundance is determined by the thermal generation and the non-thermal generation from the stau 
NLSP decay: 

^Fy = ^F,,tuh' + ^A,nontl.h'. (12) 
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First, the freeze-out density of the rehc stau can be calculated by taking the leading s— wave stau annihilation to 
photons whose thermal-averaged cross section is given by cr ~ 21:0? jrr?- with the fine structure constant a. Then, 
following the standard procedure presented in [4l[ , the non-thermal density of F\ can be estimated as 

f^i?i,nonth^' - 0-1 (400Gev) (400Gev) ' ^^^^ 

By using the analytical calculation in Ref. [1^, the density parameter of the thermal component of F\ is estimated 
to be 

Here Xq = mq/TRH, with rriq the squark mass and Trh the reheat temperature after the lightest flaton F\ decay, 

— (T^)"(T54v)lfT)"'' <-) 

with Ba the branching ratio of F\ decaying into axions. Trh should be sufficiently high for daughter axions produced 
by the decay of F\ not to induce "^He overproduction ( i?Q < 0.1 ^2^) and must be also larger than mf/2b to thermally 
produce the NLSP staus. The quantity e^^^^ denotes the squark-quark-flaxino mixing which becomes the largest for 
the stop having an order-one top Yukawa coupling hf. Sgqp-^ — hfNo4 with iVo4 ~ SiCp ~ 0(v/Fa). 



V. THE FLATON AND FLATINO MASS SPECTRUM 



By assuming the superpotential given in we have a scalar potential including soft supersymmetry breaking 
terms, 

^.oft = f4fp^Q + ^-c- (16) 

Mp 

Then the flaton masses are given by 

mp,, = -^\/3(12 - + 2^^(12 + e) ± |12 - C\Vx^ + 42a;2 + 9, (17) 
v2 

mp, = (18) 
where jl = fif /h and ^ = —Af /{ffio) with /io = (P) (Q) /2Mp Then, the flatino masses are represented by 

mp^ ^ = ( + 1 ± 1) . (I9) 

Allowed parameter regions for x and ^ of our interest (i.e, parameters where the flatons cannot decay into flatinos, 
"^-Fi,2,-F' < 2m^^ ^) are given in (28) in Ref. [11]. 



VI. RESULTS 



We present our results on the DM abundance and the stau lifetime by taking some typical parameter sets of the 
flat DFSZ model. Following the notations and the parameters presented in the previous sections and Ref. [1^, we 
take a set of X = 4, ^ = 13, f/h = 1/24, A^/fJ, = -1, Af/^i = 13/24, tan/3 = 3, and M2 = 2Mi. In this paper we 
show two cases of the mass parameters 1) ji — 287 GeV {ji — 12) and Mi ~ 130 GeV, and 2) ji = 813 GeV {jl = 34) 
and Ml = 530 GeV. 

We can calculate the mass spectrum in GeV unit of the flatons F1.2, the pseudoscalar flaton F' , and the flatinos 
Fi_2 which is given as 

r n r [216,175,162,184,112] (A = 12) , 

Lmf',mf,,mf,,m^^,m^j - I jg^g gQQ^^g2^525,320] (A = 34) . 
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FIG. 1; Allowed regions of flaxino dark matter. The axes correspond to the masses of the stau and the squark (stop). The 
other parameters are as described in the text, with the choices /i=12, x—4, ^=13, tan/3=3, and Fa = 10^" GeV. The region 
between the two solid lines satisfies the WMAP constraints on Qlsp given in Eq. ([l}. The lightest flatino LSP and the stau 
NLSP can be realized in the region between the two dot-dashed lines. The intersection of these regions is allowed. 



In each case, by using Eq.ljTS]) the reheat temperature after the thermal inflation (from the decay of the lightest 
flaton Fi) is given by Trh = 39 ( Trh = 189 GeV) for fi ^ 12 {fi = 34). 

In addition, we can diagonalize the 6x6 mass matrix by using the approximated 5x5 method discussed in Sec. Hill 
and get the mass spectrum of four neutralinos plus two flatinos. Then we can calculate the mass of the lightest 
neutralino xo (or A^i) and get 1) m^^ — 120 GeV and 2) 525 GeV, respectively. 

In Figs. [T] and [2] we plot the observationally-allowed region which is enclosed by two solid lines in (mf,TOq) plane 
for /i=12 and 34, respectively. Here we have Fa = 10^° GeV. In both figures, the vertical area enclosed by solid lines 
comes from the non-thermal contribution and the narrow horizontal area from the thermal contribution which is very 
sensitive function of Xq. The vertical band between two dot-dashed lines denotes a region where the lightest flatino 
(flaxino) LSP and the stau NLSP are realized. In Fig. [2l we see that the flaxino dark matter naturally agrees with 
observations in broad parameter spaces of rrig at the weak scale. 

In Figs. [3] and 21 we plot the lifetime of stau as a function of stau mass in case oi p, = 12, and 34, respectively. The 
thick (thin) solid line denotes the case of the mixing angle Sf = 1(0). Here we have taken Fa = 10^° GeV to minimize 
the decay length. In case of /2 = 12, the decay length is shown to be around 1 m and thus the stau NLSP decay to the 
flaxino LSP will be possibly searched for at LHC. Note that the lifetime scales (x F^. From these figures we see that 
the lifetime is sufficiently shorter than lO'^ sec excluding the special case of degenerate masses mf — mp_^^ ttIt (~ 1.8 
GeV) . Hence this scenario is not constrained by the ^Li overproduction through BBN if we adopt the standard range 
of Fa = 10^° GeV - 10^^ GeV. In other words, even lager value Fa ~ lO^'' GeV can be allowed from the viewpoint 
of BBN. Such a high value might not lead to the overclosure of the universe by the axions because of a sufficiently 
low reheating temperature and large entropy productions after the thermal inflation [2^ . See also Ref . [i^l for recent 
detailed analyses of the thermally-produced axions for a relatively low reheating temperature. 



VII. CONCLUSIONS 



The lightest supersymmetric particle (LSP) is a dark matter candidate. It is usually said that the LSP will be either 
the lightest superpartner of a Standard Model particle (so called LOSP), the gravitino or the axino. For cosmological 
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FIG. 2: Same as Fig. [T] but for fl=3A. 




FIG. 3: Stau lifetime decaying into the lightest flatino. The thick (thin) solid line denotes the case of the mixing angle Sf — 1 
(sf = 0). Here we have adopted Fa = lO^^GeV and the electroweak scale v — 264 GeV. Note that m^^ =112 GeV, and the 

mass of the lightest neutralino A''i (m^p = = 120 GeV) is smaller than that of F2 {rrip =184 GeV) in this case. 
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FIG. 4; Same as FigHbut for /i=34 {^i = 813 GeV) and Mi = 530 GeV. Then, mp = 320 GeV, and mf, ~mp = 525 GeV. 



constraints and collider signals of the LSP, it matters what is the NLSP. In this paper we revisited the axino-like LSP 
arising in the flat DFSZ axion models for the case of the stau NLSP. 

Following [28| , we assumed that the spontaneous breaking of Peccei-Quinn symmetry comes only from supersym- 
metry breaking. Then all of the PQ particles except the axion will have masses of order the gravitino mass. The 
scalar particles are called flatons and their superpartners are called flatinos. 

In this scenario the saxion is a linear combination of flatons and the axino a linear combination of flatinos, but they 
are not mass eigenstates and have no special status. So the LSP candidate from the PQ sector is not the axino, but 
the lightest flatino which we have dubbed the flaxino. We have explored the simplest version of this scenario, which 
has two flatinos and assumes gravity mediated supersymmetry breaking. 

We have found that the stau NLSP decay to the flaxino LSP is fast enough to maintain the standard predictions 
of Big Bang Nucleosynthesis. The big difference from the gravitino dark matter scenario is that the flaxino/axino 
coupling (~ ^/Fa) is much larger than the gravitino coupling (~ 1/Afp), which makes the stau NLSP decay much 
faster. In the DFSZ axion model, the stau-flaxino/axino coupling appears at tree level as a consequence of the /i term 
mixing between the higgsino and flaxino/axino. This makes the NLSP decay more efficient compared with the KSVZ 
axion models. 

In a specific DFSZ model whose parameters are constrained by various cosmological considerations, we identified 
the region of the stop and stau masses which is compatible with the required dark matter density. Within the 
conventional bound on the axion scale Fa ^ 10^^ GeV, the stau decay to the fiaxino LSP never becomes dangerous 
for BBN. Furthermore, with the chosen parameter sets in the first case (/i — 12), the stau decay length turns out to 
be 0(1) m so that a fair number of stau decays can be captured in colliders when the axion scale is close to its lower 
bound. Fa ~ 10^" GeV. In the second case (/2 — 34), the flaxino LSP is naturally realized in the broad parameter 
space of the squark mass. Although the decay length is longer than 0(10) m, the additional detectors proposed by 
[11, [3^ will catch the staus even in this case. 

Let us finally remark that, in our scheme, the unstable gravitino decaying to ordinary superparticle or flaxino does 
not contradict with BBN as the primordial gravitino abundance will be diluted away by the entropy production after 
thermal inflation caused by the flat direction in the model. A gravitino even lighter than the fiaxino is also allowed 
since the flaxino decay to the gravitino and the axion causes no trouble with BBN. 
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